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ABSTRACT 

We present results from a 244 ks NuSTAR observation of 3C 273 obtained during a cross-calibration 
campaign with the Chandra^ INTEGRAL^ Suzaku^ Swifts and XMM-Newton observatories. We 
show that the spectrum, when fit with a power-law model using data from all observatories ex¬ 
cept INTEGRAL over the 1-78 keV band, leaves significant residuals in the NuSTAR data between 
30-78 keV. The NuSTAR 3-78 keV spectrum is well-described by an exponentially cutoff power-law 
(r = 1.646±0.006, Ecutoff = 202^11 keV) with a weak reflection component from cold, dense material. 

There is also evidence for a weak {EW = 23± 11 eV) neutral iron line. We interpret these features as 
arising from coronal emission plus reflection off an accretion disk or distant material. Beyond 80 keV 
INTEGRAL data show clear excess flux relative to an extrapolation of the AGN model fit to NuSTAR. 

This high-energy power-law is consistent with the presence of a beamed jet, which begins to dominate 
over emission from the inner accretion flow at 30-40 keV. Modeling the jet locally (in the NuSTAR 
+ INTEGRAL band) as a power-law, we find the coronal component is fit by Tagn = 1.638 ± 0.045, 

Ecutoff = 47 ± 15keV, and jet photon index by Tjet = 1-05 ± 0.4. We also consider Fermi/LAT ob¬ 
servations of 3C 273 and here the broad-band spectrum of the jet can be described by a log-parabolic 
model, peaking at ~ 2MeV. Finally, we investigate the spectral variability in the NuSTAR band and 
find an inverse correlation between flux and T. 

Subject headings: quasarsdndividual (3C 273) - X-rays: individual (3C 273) 


1. INTRODUCTION 


At a redshift oi z = 0.158 (Schmidt 


1963), 3C273 is 


the nearest high luminosity quasar and has been exten¬ 
sively studied at all wa velengths since its discovery in 


1963 (for a review, see Courvoisier 1998). It is radio- 


loud and hi ghly variable across nearly all energies (Soldi 


et al.l 2008), with a jet showing apparent superluminal 


motion. 

At radio to millimeter and at y-ray energies. 


flares 


bump) which is spectrally complicated and suggestive of 
two independently varying components: a more rapidly 
varying component originating from thermal reprocess¬ 
ing from an accretion disk, and a slower component that 
cou ld be synchrotron emission unrelated to the radio-mm 
jet ( Paltani, Courvoisier fc Walterjl998' ). As observed in 


many other ACN, there is a soft-excess in the low-energy 
(< 2 keV ) X-ray band pos sibly due to Comptonized UV 
photons (Page et al.|2004[). Correlations between the UV 


and low-energy X-rays (Walter &: Courvoisier|1992 Page 


et al.||20d4 ) have been noted in the past, supporting this 
interpretation. However, recent observations have failed 
to detect correlat ed variability (Chernyakova et al.|2007 


Soldi et T]|2008), leaving the interpretation of the opti¬ 


cal/UV excess uncertain. 

There is evidence of an intermittently weak iron line 
in the X-ray spectrum, which appears to be bro ad (a 
0.6 keV, EW ^ 20-60 eV), occasionally neutral (Tur^r[ 


et al. 1990 Page et al. 
and sometimes 


2004 Grandi & Palumbo 20041 


(iTiirler, Courvoisier & Paltani 

2000 

Abdo et al. 

2010 

1 . 


Kataoka et al. 

2002 

trace its varia 

bility 


ionized (|Yaqoob & Serlemitsos| |2000 
The line is too faint to reliably 


instrumentation. 

Above 2keV and up to MeV energies, previous obser¬ 
vations report a hard power-law spectrum, as is common 
for jet-dominated AGN (blazars). Over the 30 years that 
3C 273 has been reliably monitored, there appears to be 
a long term spectral evolution underlying the short term 
variations. The source was in its softest observed state 
in June 2003 (photon index P ~ 1.82 ± 0.01), a value of 
AP ~ 0.3-0.4 above what was measured in the 1980’s 
(P ~ 1.5), and since then the source has hardened again 
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to a value of F ~ 1.6 — 1.7 ( Chernyakova et ^|2007 ). 

Finally, 3C 273 is a strong 7 -ray emitter; it has been 
routinely de tected by Fermi, showing several interesting 
7 -ray flares (Abdo et al.|2010 ). Broad-band observations 
including F ermi and IJy 'lEURAL reported by |Esposito| 
et al. (2015) imply that the broad-band spectrum of the 


jet emission is not a single power-law, but instead, it can 
be described by a log-parabola model peaking around 
3-8 MeV. 

Because of the presence of the soft excess and the iron 
line, it has been postulated that the X-ray emission arises 
from a mix of thermal and non-thermal processes: an 
AGN component from the inner accretion flow (i.e. emis¬ 
sion from the corona possibly accompanied by disk reflec¬ 
tion) and a jet component. Grandi & Palumbo (2004) 
fitted the broadband spectrum from BeppoSAX with a 
cold reflector irradiated by an isotropic coronal X-ray 


sourc e for the AGN component (Magdziarz & Zdziarski 
19951 and a power-law for the jet component, 'i'hey were 


able to obtain reasonable fits to a restricted model where 
the coronal continuum photon index and cutoff energy 
were held fixed. A search fo r the disk reflection hump 
at ~30keV was presented in Chernyakova et al. (2007) 
using XMM-Newton and INTEUKAL data. I'hey esti- 
mated that about 20 % reflection could be allowed, but 
no clear results could be obtained because of the non¬ 
overlapping energy bands between the instruments and 
the gap from 10-20 keV, making cross-calibration uncer¬ 
tain. 

The two compone nt scenario i s supp orted by variabil¬ 
ity studies made by Soldi et al. (2008), who found that 


the amplitude of the variations increases steeply above 
20 keV with no convincing correlation between the vari¬ 
ations above and below 20keV. As of yet no clear sig¬ 
nature above lOkeV of a disk related component in the 
form of a reflection feature and/or up-scattering of the 
thermal disk photons by a corona has been found. We 
show that this is no longer the case and present here 
overlapping contemporaneous observations of 3C 273 ob¬ 
tained by the six observatories: Chandra, INTEGRAL, 
NuSTAR, Suzaku, Swift, and XMM-Newton. We place 
particular emphasis on the NuSTAR observation, which 
was about 6 times longer than the rest and for the first 
time allows a study of the 3-78 keV region without cross¬ 
calibration concerns. 

2 . DATA ANALYSIS 


iWeisskopf et al. 

2002), INTEGRAL ( 

Winkler et al. 

20U3), NuSTAR ( 

Harrison et al. 20131) 

^uzaku (Mit- 

suda et al.|2007D, Swift (IGehrels et ai.|2UU4 

), and XMW 

Newton 1 

Jansen et al. 

2001|) observed 3G273 as part of 


tional Astronomical Gonsortium for High Energy Cali¬ 
bration (IACHECQ. Figurej^shows the overlapping GTI 
times and length m each observation and Table[^lists the 
observation ID’s and exposure times for the respective 
observatories. 

For Chandra we used CIAO 4.6.1 and CALDB 
4. 6 .1.1. The data were taken with gratings configu¬ 
ration ACIS-I-HETG and reprocessed using the CIAO 
chcLndra_repro reprocessing script. We combined orders 

^ http://web.mit.edu/iachec/ 


1-3 for the HEG and MEG arm separately, and binned 
the data at 30 counts. 

For INTEGRAL we used the standard OSA 10.0 data 
reduction package. We only used data from IBIS/ISGRI 
(Ubertini et al.||2003 Lebrun et al. 2003) as data from 
JeM-X and SFI did not allow to constrain the spectral 
shape significantly. We extracted the IBIS data on a sci¬ 
ence window (ScW) by ScW basis using our own scripts, 
which are based on the IBIS Analysis User Manual pro¬ 
cedure. We combined the spectra using spe_pick ac¬ 
cording to observation date. 

For NuSTAR we used HEAsoft 6.15.1 and CALDB 
20131223. The data were processed with all standard 
settings and source counts extracted from a 30" radius 
circular region. Background was taken from the same 

For Suzaku we used CALDB: HXD (20110913), XIS 
(20140203) and XRT (20110630). The observation was 
taken in l/4 window mode, and we used 100 arcsec radius 
circular regions for the FI detectors (XIS0,3) and a 140 
arcsec region for the BI detector (XISl), such that the 
regions were centered on the source, but were restricted 
to the operational portions of the detectors. 

For Swift we used HEAsoft 6.15.1 and XRT CALDB 
2014-02-04. The data were taken in ‘PHOTON’ mode 
and were reduced using xrtpipeline. Spectra were ex¬ 
tracted from an annulus region, inner radius 5" and outer 
radius 30" to correct for pileup. The two observations 
were combined and the spectra binned at 50 counts. 

For XMM-Newton we used SAS v. 13.5.0 with CALDB 
2014-01-31. The data were taken in ‘Small Window’ 
mode, and to corrected for pileup in the MOS we ex¬ 
cised counts form an annulus region with inner radius of 
15" and outer radius of 45". For the PN we extracted 
from a circular region of radius 45". 

We analyzed the high-energy 7 -ray data from the 
Fermi Large Area Telescope (LAT), using the software 
package ScienceTools vlOrOpS and the latest calibration 
standard Pass 8 . Since 3C 273 was found in a rela¬ 
tively low 7 -ray state, this required integrating the 7 - 
ray flux over a timescale of 50 days (MJD 56100-56150), 
which is significantly longer than the 6 -day NuSTAR 
campaign. We used the instrument response function 
P8R2_SOURCE_V6 (front and back), including galac¬ 
tic diffuse emission model glLiem_v06, isotropic back¬ 
ground model iso_P8R2_SOURCE_V6_v06, and back¬ 
ground point sources within 15 deg from 3C 273 taken 
the 2FGL catalog. Events were extracted from a region 
of interest within 10 deg from 3C273, we applied stan¬ 
dard selection cuts for the SAA avoidance, and the zenith 
angle cut z < 100 deg. 

3. SPECTRAL FITTING 

For all the data anal ysis we use th e XSPEC version 
12.8.2 analysis software ( Arnaud|1996 ), and present 90% 
confidence limits unless otherwise stated. 

3.1. Broad-band Fitting 

To investigate the overall broadband spectrum we fit 
data to all observatories except INTEGRAL (to be dis¬ 
cussed in the next section) for an interval with reasonable 
temporal overlap. Figure shows the observation times 
for the various missions. \Ve chose the interval span¬ 
ning the Suzaku observation for extraction of data since 
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TABLE 1 
Observation Log 


Instrument 

Start Time 

OBSID 

Exposure (ks) 

Chandra 

2012 Jul 16 UT 11:04:29 

14455 

30.0 

NuSTAR 

2012 Jul 14 UT 00:06:07 

10002020001 

244.0 

Swift 

2012 Jul 16 UT 10:24:59 

00050900019 

13.0 

Swift 

2012 Jul 17 UT 00:50:59 

00050900020 

6.9 

Suzaku 

2012 Jul 16 UT 08:08:54 

107013010 

39.8 

XMM-Newton 

2012 Jul 16 UT 11:59:23 

0414191001 

38.9 

INTEGRAL 

2011 Dec 05 UT 02:27:04 

111^1128 (321 ScW) 

689.5 

INTEGRAL 

2012 Jun 08 UT 04:08:54 

1178, 1180-1183, 1188 (66 ScW) 

131.6 

INTEGRAL 

2012 Jul 14 UT 22:18:59 

1191-1192 (99 ScW) 

68.4 

Instrument 

Start Time 

Stop Time 

Exposure (days) 

Fermi 

2012 Jun 22 UT 00:00:00 

2012 Aug 11 UT 00:00:00 

50 


‘‘Revolutions. 


integral ^ ^ ' ' ' ' " 

Chandra hh hb 

NuSTAR ■■■■ ■■ ■■ 

Suzaku 

Swift 

XMM-Newton 

IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIN 

lllllllllllllllllll 

lllllllllllllll 


■I 


-2 0 2 
Time (Days since UT 00:00:00 2012 July 17 ) 


Fig. 1.— Length and coverage of the six overlapping 3C 273 ob¬ 
servations. We do not enforce strict simultaneity between the low 
earth orbit observatories NuSTAR, Suzaku, Swift, and extract the 
full range for all observatories - except for INTEGRAL - between 
the Suzaku start and stop times as indicated by the shaded region 
in Figure]^ 


there is reasonable overlap of part of this interval with 
all observatories except INTEGRAL. Due to the relative 
phasing of the South Atlantic Anomaly (SAA) passages 
and occultation periods among the low Earth orbit ob¬ 
servatories {NuSTAR, Suzaku, and Swift), we do not en¬ 
force strict simultaneity within this window. The NuS¬ 
TAR lightcurve in Figure (binned at 25 ks, live-time 
and vignetting corrected) snows that during the Suzaku 
window, marked by the shaded region, the flux changed 
only gradually by ^10%. We therefore do not expect 
sudden flux changes to have occurred between occulta- 
tion/SAA periods. We extracted the full time range for 
all instruments except NuSTAR where we truncated the 
observation by applying the Suzaku start and stop times. 
Because the majority of the INTEGRAL observation lies 
outside the near simultaneous window, where the flux is 
significantly different, we excluded INTEGRAL from this 
broadband fit. 

We fitted an absorbed power-law niodel using Wilms 
abundances ( Wilms, Allen fc McC ray||2000 1 and Verner 
cross-section^ Verner et al. 199b| , freezing the hydro- 
gen column to the Galac tic value of 1.79x10^° cm“^ 
(Dickey & Lockman 1990). We keep this value fixed in 
the remainder ot the hts. To avoid complications with 
the soft excess, we limited our fitting range between 1- 
78keV, though we note that during these observations 
the soft excess appears to have been very modest. We 
allowed all normalizations to float; the relative cross¬ 
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Fig. 2.— Light-curve in different energy bands for NuSTAR. Bin 
width is 25ks and the light-curve has been live-time and vignetting 
corrected. The shaded region marks the start and stop of the 
Suzaku observation. 
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normalization terms and the relative slope errors between 
instruments are discussed in detail in Madsen (in prep) 
by the lACHEC consortium. 

Figure shows the results of the best fit power-law: 
r = 1.647 ± 0.003. Formally the fit is good, with 
Xrod = 1-018 (7815 dof). However, visual inspection of 
the high energy residuals above 20 keV in NuSTAR shows 
obvious systematic deviations from a power-law. This 
good formal fit is due to the very large number of bins 
between 2-8 keV. Replacing the power-law with a cutoff 
power-law, cutoffpwrlw (XSPEC), reduces the residu¬ 
als, and we find T = 1.624 ± 0.006, Acutoff = 29ll“ keV 
(Xred = 1.013). To examine any further features we fit¬ 
ted a cutoff power-law to NuSTAR alone, ignoring the 
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energy ranges 5-8 keV and 10-50 keV and plot the NuS- 
TAR ratio of the model to data alone in Figure]^ The 
shape of these residuals is reminiscent of reflection of the 
primary continuum off dense material, with an iron line 
and a Compton hump. The NuSTAR calibration uncer¬ 
tainties are ± 1% up to lOkeV, ±2 % up to 40keV, and 
5-10% above (Madsen et al. 20151. While the excess is 
at a level only somewhat greater than the calibration un¬ 
certainties, the systematic shape is not characteristic of 
calibration errors, and we associate it with a very weak 
reflection component. 

We then fitted all the instruments with the combina¬ 
tion of a primary cutoff power-law plus reflection from a 
plane parallel sl ab of neutral material (the pex rav model 
in XSPEC; see (Magdziarz & Zdziarski 1995)). The jet 
angle and hence disk orientation h as been estimated from 
radio observations to be ~ 10° (Abraham & Romero 


1999), though fits to the UV-conti nuum and iron line 


su ggest it could be as high as ~ 60° (Yaqoob & Serlemit- 


sos 


2000). Assuming the material to be associated with 


the accretion disk, we tested the fits for both angles, but 
since we found no appreciable effect on the fit, we fixed 
the value at 35°. We used solar abundances and the best 
fit is achieved for T = 1.646±0.006, Ecutoff = 202+34 keV, 
with a weak relative reflection component with ratio of 
direct to reflected flux oi R = 0.15±0.05 (x^ed = 1-006). 
We show the ratio of model to data of this fit in Figure 
[H A Monte-Carlo simulation shows the likelihood of a 
reflection of this order to be in excess of 99.9%. 


3.2. Fitting NuSTAR and INTEGRAL 

We investigated joint NuSTAR and INTEGRAL ob¬ 
servations over a broader time range, dispensing with 
the low energy-instruments. We applied the INTEGRAL 
GTI’s from the July observation to the NuSTAR data 
set. Fitting the data with a pexrav alone results in IN¬ 
TEGRAL diverging above 80keV. A power-law provides 
a reasonable broad band fit, but with strong NuSTAR 
residuals between 20-78 keV as shown in Figure Fol¬ 
lowing the idea that the X-ray spectrum is in reality a 
superposition of an AGN component (i.e. coronal emis¬ 
sion plus disk reflection) with a jet component, we fitted 
the data with the pexrav model for the AGN compo¬ 
nent and a flux pegged power-law, pegpwrlw (XSPEG), 
for the jet. 

The two models prove to be degenerate. To partially 
break this degeneracy we investigated three INTEGRAL 
data sets taken in December 2011, June 2012, and July 
2012. We fitted the spectra with pegpwrlw from 30- 
250 keV. The individual fits are summarized in Table [2] 
with the spectral energy distributions shown in Figure 
The spectral slope changes with time, but the flux 
difference decreases with increasing energy and converges 
at ~ 120 keV. We also draw attention to the dip in the 
June and July observations between 60-80 keV; the error 
bars are large, but this could be interpreted as the AGN 
turn-over seen by NuSTAR. 

Based on the fact that for different slopes of the spec¬ 
trum in INTEGRAL the flux between 80-150 keV re¬ 
mained almost constant, we decided to freeze the flux 
of the jet component between 80-150 keV to 92 x 10“^^ 
erg cm“^ s“^ as obtained from the July observation. Ad¬ 
ditionally, we freeze the relative reflection to R=0.15 to 
further limit the degeneracy, but allowed all other pa¬ 


TABLE 2 

INTEGRAL Spectral fits: pegpwrlw 


Name 

r 

Flu> “ 


Flm 



X?p,i/dof 

December 2011 

1.76 ±0.09 

118 ±5 

91 ± 7 


4.90/6 

June 2012 

1.88 ±0.16 

128 ± 10 

94 ± 14 

4.66/6 

July 2012 

1.53 ±0.25 

105 ± 13 

92 ±20 

1.49/6 


‘‘Energy range 40-100 keV in units of 10^^ erg cm^ s 
'‘Energy range 80—150keV in units of lO'^^ erg cm'^ 


rameters to vary. 

Table lists the best fit parameters. We explore the 
degeneracies between Fpexravi and rp 0 gp^i.i^ in 

Figure The spectral energy distribution of the best 
fit with the two model components is shown in panel A. 
Panel B shows the flux ratio between the AGN/jet for 
three different 80-150 keV jet fluxes (Fjet = 72, 92, and 
112x10“^^ erg cm“^ s“^). These fluxes span the ±lcr 
level errors of the INTEGRAL July fit (Table [^, and 
for this range of fluxes the point where the jet crosses 
over the AGN component occurs between 30-50 keV. At 
6 keV the AGN component is about an order of magni¬ 
tude above the jet. Contours of Fpexrav versus Ecutoff 
are shown in panel C, and Fpexrav versus Fpegpwriw in 
panel D. While these parameters are clearly coupled, 
they bound a well defined space, predicting the AGN 
photon index to be between 1 . 6 - 1 .7 with a cutoff energy 
of 30-70 keV. The jet photon index lies somewhere be¬ 
tween 0.5-1.5. 

To get a better constraint on the slope of the jet, we 
extracted a high-energy y-ray (> 100 MeV) spectrum of 
3G 273 from the Fermi/LAT data. We assumed that the 
broad-band SED of the jet component is a log-parabola 
N{E) = 7VLogp(E/Fo)[-“-'^^°s(i5/i5o)]^ which has the ad¬ 
vantage of adding only one parameter to the overall 
model. Parameter a corresponds to the local photon 
index at the photon energy of Eq, and here we choose 
Eq = 100 keV. We added this model to the best-fit 
pexrav model from Table 1^ which we kept fixed, and 
we fitted it jointly to theNuSTAR, INTEGRAL and 
Fermi/LAT data (Figure]^. We obtained a good fit 
with Xred = 1-117, with the results a = 1.428 ± 0.011 
and /? = 0.091 ± 0.002. This places the SED peak of the 
jet component at Euea> = (2.28 ± 0.21) MeV, consistent 
with the results of Esposito et al. (2015). 

We note that the INTEGRAL data points are system¬ 
atically above the fit. But a perfect fit is not expected, 
since 3G273 was found in a relatively low y-ray state, 
which required integrating the y-ray flux over a timescale 
of 50 days, far longer than the 6 -day NuSTAR campaign. 
Short time scale variability might therefore have put the 
jet in a higher flux state during the NuSTAR period than 
on average. The fit, however, shows that the inferred 
slope of the jet is indeed reasonable. 


3.3. The Iron Line 

Although the detection of an iron line in 3C 273 has 
in the past been firm, it remains faint and at times even 
absent. To maximize our signal we use all of the 244 ks of 
the NuSTAR observation and just NuSTAR alone. We 
restricted the energy range between 3-10 keV. In this lim¬ 
ited band an absorbed power-law provides a good fit, but 
with a sli^t excess between 5-7 keV as can be observed 
in Figured We added a gaussian component with a fixed 
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Fig. 3.— Top panel: Simultaneous broadband absorbed power- 
law fit between 1—78keV to all observatories except INTEGRAL. 
Cross-normalization constants are allowed to float and the column 
is fixed to the Galactic value of Wh=1-79x 10^*^ cm~^. The best 
fit photon index is F = 1.647 di 0.003. Bottom panels: Ratio 
of a power-law (middle panel) and the pexrav model (bottom) 
shown for NuSTAR alone, pexrav fit: F = 1.646 di 0.006, EcutofF = 
2021 * 13 ^ keV, and relative reflection oi R = 0.15 di 0.05. 



Energy (keV) 


Fig. 4. — Ratio of NuSTAR data alone to a cutoff power-law 
model ignoring the energy ranges 5—8 keV and 10-50 keV. The NuS¬ 
TAR calibration errors are at the ±1% level, which accounts for 
the spectrum’s scatter below 5 keV. The iron line and Compton 
reflection hump are very modest but systematically different from 
the known calibration errors. 

rest energy of 6.4 keV and the quality of fit improves from 
X^/dof = 419/345 without to x^jdoj = 398/343 with. 
A Monte-Carlo simulation shows this feature to be sig¬ 
nificant in excess of 99.97%. 


The width of the line is broad, CTwidth = 0.65 ± 0.3 keV, 
and we measure an equivalent width of EW = 23± 11 eV, 
which is w ithin the envelope of what has been previously 


reported (Yaqoob & Serlemitsos 20001 Kataoka et al. 

2002 

Page et al.||20U4 

Crandi & l’alumbo||2004 ). 


We note that the magnitude of the iron emission is 
close to the calibration limit (±1%), and as a separate 
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Fig. 5. — Residuals of the data to model of the fits between NuS¬ 
TAR and INTEGRAL for a pexrav (top) and pegpwrlw (bottom) 
fit. Only NuSTAR FPMA is shown for ease of viewing. 
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1 ® 
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Fig. 6. — Power-law fits to INTEGRAL data. Flux differences 
decrease with increasing energy and appear to converge around 
~ 120 keV. Fit parameters are recorded in Table 

test we fixed the line energy to 6.4 keV and allowed 
the redshift to vary. We recovered the correct redshift 
within errors. We also replaced the line with a relativis¬ 
tic diskline model but did not find significant improve¬ 
ment to the fit. 


3.4. Continuum Models 

As discussed in the previous sections, the 3C 273 con¬ 
tinuum data can be well fit with an absorbed cutoff 
power-law with weak reflection (a pexrav model) up to 
the NuSTAR upper energy limit at 78keV. However, in¬ 
clusion of the INTEGRAL data reveals that the best 
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Fig. 7. — INTEGRAL and NuSTAR fit. A) fit from Table ( black/red = FPMA/FPMB, green=ISGRI). B) Ratio between jet and AGN 
component for 3 flux settings of the jet (given in units of 10"^ erg cm~^ C) Contours of input spectrum of the pexrav. D) Contours 

of the photon index of the AGN against the jet component. 


TABLE 3 

INTEGRAL/NuSTAR June: tbabs*const*(pexrav+pegpwrlw) 


pexrav (AGN) 


Xref 

0.973 for 1478 dof 

TylGiV 

1.638 ±0.045 

Cutoff Energy 

47 ±15 keV 

Relative reflection 

0.15 (frozen) 


pegpwrlw (jet) 


F jet 

Flu>p| (80-150 keV) 

1.05 ± 0.4 

92 (frozen) 

constant FPMA 

1.0 (frozen) 

constant FPMB 

1.04 ±0.01 

constant ISGRI 

0.95 ±0.11 

1 9 —9 —1 

erg cm ^ s ^ 



Fig. 8.— Broad-band SED for 3C 273 where the AGN component 
(pexrav) has been kept fixed to the values from Tablel^and the jet 
fitted with a log-parabola that peaks at Epeak = 2.28i 0.21 MeV. 


I 


fit between the two observatories requires an additional 
power-law for the jet component. Unfortunately, INTE¬ 
GRAL does not cover the entire NuSTAR observation, 
and the rise of the jet component is not visible in the 
NuSTAR data alone. We therefore proceed with our 
investigation of the continuum by using the full 244 ks 
NuSTAR observation, but make the assumption that the 
jet component is constant for the entire time range and 
takes on the values given in Table We realize this to 
be an approximation and therefore to contrast we also 
fit the the data without the jet component to bracket 
parameters. 

The results of the continuum fits are shown in Table 
^ and the residual of data to model in Figure We 
rst fit the continuum with a pexrav model with and 
without a gaussian iron line inc luded and compare with 
a pexmon (Nandra et al. 2007), which self-consistently 
includes the iron line, i'he two models yield very similar 
parameters with a photon index T = 1.63 ± 0.01 and 
Ti^cutoff ~ 260±35keV for the model without the jet, and 
r = 1.66 ± 0.01 and Ffcutoff ~ 50 ± 2keV with the jet. 
The relative reflection for both cases is around R ~ 0.1. 

To probe a more physical model for the coronal con- 
tinuum, we fitted with a CompTT model ( Titarchuk|1994 1 
and included the reflection component using pexmon. 
The CompTT model calculates the emergent spectrum re¬ 
sulting from Comptonization in an electron plasma, tak¬ 
ing into account relativistic effects. It has models for a 
disk and spherical plasma geometry, and we adopted the 
disk geometry. We assumed solar abundances, fixed the 
inclination angle at 35°, and used as soft photon tern- 
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ponent to the values found for the pexmon fit, check¬ 
ing that the shape of the direct spectrum of the pexmon 
model matched that of the CompTT direct spectrum, and 
included only the reflection part of the model. 

Table shows the results from fitting the CompTT disk 
model to the NuSTAR data with and without a jet com¬ 
ponent (with jet parameters again fixed to match the 
INTEGRAL E > 80 keV data as provided in Table [^. 
The inclusion of the jet component significantly changes 
the best fit parameters. Without the jet, we find val¬ 


ues of the plasma temperature of fcTe = 247^g|keV, 
while with the jet included the temperature decreases to 
/cTe = 12.0 ± 0.3 ± 1 keV and increases the optical depth 
from Tp = 0.15lg;g4 without jet to Tp = 2.77 ± 0.06 with 
the jet. As shown (Figure the ht without jet is ex¬ 
cellent and shows no indication of a jet component. The 
fit with jet displays residuals at high energies, indicating 
that our assumption of a constant jet is not correct, or 
that the constraints from INTEGRAL are not sufficient 
to properly constrain this component. Allowing the jet 
to fit freely we can find better solutions, but these results 
push the CompTT model up against its valid phase space 
and are thus questionable. We can, however, place an 
upper limit of kT^ ^ keV on the electron tem¬ 

perature, since inclusion of the jet in any form decreases 
the value from the no-jet scenario. Conversely the op¬ 
tical depth has a lower limit of Tp = 0.33 ± 0.1. Given 
that the uncertainties in the CompTT physical parameters 
are completely driven by the uncertainties in the level 
of the jet flux, we do not investigate alternative coronal 
geometries further. 


We attempted to fit with a CompPS model (Poutanen 
& Svensson|1996|) , which is a coronal plasma rnodel that 
includes reflection off a cold disk in the manner of the 
pexrav, but the returned fits were within a regime of 
high optical depth and low temperature where the model 
is unfortunately not c onsidered reliable due to the large 
amount of scattering (jPoutanen &: Svensson||1996). 


3.5. Variability 


The NuSTAR lightcurve (Figure]^ shows that 3C273 
went through flux changes of 10 -30% during the ob¬ 
servation. We calculated the noise subtracted variance 



We investigated the variability of the spectrum by 
decomposing the observation into 5 time bins, GTIA- 
GTIE (start and stop times are listed in Table [^, as 
illustrated in Figure [TT| The motivation behind this par¬ 
titioning is to cover the sections where the flux is roughly 
constant. We chose to continue with the pexrav model, 
and fitted both with and without the jet component, as¬ 
suming it was frozen to the values given in Table The 
two models and the parameters are summarized in Table 

El 

We find that by enforcing a jet component we obtain 
lower Fcutoff energy since more bending is required when 
the spectrum is diluted by the jet, and the jet compo¬ 
nent systematically softens the index. The jet compo¬ 
nent also systematically increases the relative reflection 
for the same reason. 

Because of the coupling between the photon index, cut¬ 


TABLE 4 

NuSTAR CONTINUUM SPECTRAL FITS. 


Parameter 

w/o jet 1 w/ jet 

pexra\|^| 

(835 dof) 

Xfed 

r 

-^cutoff 

Relative reflection 

1.042 

1.63 ±0.01 
262 ± 34 keV 
0.07 ± 0.03 

1.039 

1.66 ± 0.01 

52 ± 2 keV 

0.2 ± 0.05 

pexrav^+zgauss 

(833 dof) 

Xred 

r 

-^cutoff 

Relative reflection 

Line Energy (fixed) 

Line width cr 

EW 

1.033 

1.63 ±0.01 
260)';3g keV 
0.07 ± 0.03 

6.4 keV 
0.62+° l° keV 
15 ± 11 eV 

1.035 

1.66 ± 0.01 

52 ± 2 keV 

0.2 ± 0.03 

6.4 keV 
0.87j:g;? keV 

15 ± 11 eV 

pexmon®' 

(835 dof) 

Xred 

r 

-^cutoff 

Relative Reflection 

1.035 

1.63 ±0.01 
257/^^ keV 
0.04 ± 0.01 

1.055 

1.64 ± 0.01 

53 ± 2 keV 

0.08 ± 0.02 

CompTT ®'+pexmon 

(840 dof) 

Xied 

Geometry 

To (fixed) 

Plasma temp, kT 

Tp 

r (fixed) 

Ecutoff (fixed) 

Relative Reflection (fixed) 

1.040 

disk 

10 eV 

nsAfg keV 

Q oo±0.11 
U.OO_o gg 

1.63 

257 keV 

0.04 

1.145 

disk 

10 eV 

12.0 ± 0.3 keV 

2.77 ± 0.06 

1.64 

53keV 

0.08 


^Solar abundances. 


off energy, and relative reflection, it is hard to determine 
the correlation between flux and slope for the AGN com¬ 
ponent. To approach this question in a model indepen¬ 
dent manner, we decomposed the spectrum into four en¬ 
ergy bands, 3-10, 10-20, 20-40, and 40-78 keV, and fit 
each section with a phenomenological power-law, calcu¬ 
lating the flux for the intervals. The results are shown 
in Figure |12| The highest flux points marked with an 
open symbm is GTI A. We fitted a linear function to all 
data points (solid line) and to GTIB-GTIE excluding 

the GTI A point (dashed line). _ 

The linear Pearson correlation (Press et al.|1992 1 favors 
in all but the lowest energy range the GTlB-GTlE (PI 
series) over GTIA-GTIE (P2 series): (3-lOkeV): P1/P2 
= -0.73/-0.86, (10-20keV): P1/P2 = -0.79/-0.64, (20- 
40keV): P1/P2 = -0.64/0.14, and (40-78keV): P1/P2 
= -0.70/-0.62. In the lowest energy band we find that 
the photon index is inversely correlated with the flux, 
becoming harder for higher fluxes. This correlation per¬ 
sists through all bands when excluding GTI A, but starts 
to break down in the 10-20 keV band when GTI A is in¬ 
cluded. Since the source during GTI A was very stable 
compared to later times, it might suggest the source was 
in a different state altogether. 


4. DISGUSSION 

From the radio to y-ray bands the SED of 3C 273 is 
dominated by non-thermal radiation from the jet, how¬ 
ever, evidence for emission from the inner accretion flow 
has been argued based on the presence of a weak iron line 
and soft excess in the X-ray, and the big blue-bump in 
the optical/UV, which are all common features of AGN 































100 

Energy (keV) 


Fig. 9.— Ratio plots of data to model for the fits in Table|^ Not shown is the pexrav+zgauss since the improvement in these residuals 
is in the iron-line region covered in Figure]^ 


TABLE 5 

NuSTAR Spectral fits: tbabs x pexrav 


Name | START | STOP | F | Ecutoff (keV) | Relative reflection | 


w/o jet 


GTIA 

2012:196:00:36:14 

2012:198:01:59:59 

1.61 ± 0.02 

213715 

0.07 ± 0.04 

1.056/1360 

GTIB 

2012:198:01:59:59 

2012:199:05:46:39 

1.67 ± 0.02 

2577“® 

0.13 ±0.07 

0.993/1036 

GTIG 

2012:199:05:46:39 

2012:200:02:36:39 

1.63 ±0.03 

265770^ 

U.UO_Q Q5 

1.025/1141 

GTID 

2012:200:02:36:39 

2012:201:06:23:19 

1.65 ±0.02 

48575 ?^ 

^•^^-0.05 

0.978/1035 

GTIE 

2012:201:06:23:19 

2012:201:23:28:43 

1.63 ±0.03 

29378^ 

-0.07 

1.047/1006 


w/ jet 


GTIA 

2012:196:00:36:14 

2012:198:01:59:59 

1.64 ± 0.02 

54 ± 4 

0.19 ±0.07 

1.051/1360 

GTIB 

2012:198:01:59:59 

2012:199:05:46:39 

1.69 ±0.03 

45 ± 5 

0.33 ±0.13 

0.996/1041 

GTIG 

2012:199:05:46:39 

2012:200:02:36:39 

1.65 ±0.04 

48 ±6 

0.20 ±0.13 

1.038/956 

GTID 

2012:200:02:36:39 

2012:201:06:23:19 

1.68 ±0.03 

55 ± 7 

0.16 ±0.11 

0.977/1035 

GTIE 

2012:201:06:23:19 

2012:201:23:28:43 

1.66 ±0.04 

55 ±8 

0.18 ±0.13 

1.039/893 


(iKataoka et al. 

20021 Grandi & Palumbo 

2004 

Pietrini 

& 'I’orricelli-Ciamponi 

2(JU8|). Evidence tor hard X-ray 


emission from tfie corona in the form of a cutoff power- 
law and reflection features above 10 keV have been elu¬ 
sive, perhaps due to the available instrument sensitivities 
or because the source was observed in states with high 
jet flux, which could have obscured the AGN signatures. 

The NuSTAR data show that the 3-78 keV spectrum 
cannot be fit simply by an absorbed power-law, but the 
spectrum requires a roll over above ^20keV. Modeling 
the NuSTAR data alone we find evidence for a weak re¬ 
flected continuum consistent with reflection from a cold 
accretion disk or distant material. We confirm the pres¬ 
ence of a weak iron line and assuming it is neutral, we 
measure its parameters to be cr = 0.65 ± 0.3 keV with an 


EW = 23±11 eV, consistent with previous measurements 
(lYaqoob fc Serlemitsos|2000l Kataoka et al. 2002 Page 


et al.|2O0l) ' Urandf^ Palumbo|2004| ) . it we include a re- 


tlection component we can account for the weak line and 
excess above 10 keV (the compton hump) with a reflec¬ 
tion fraction R= 0.07 ± 0.03 and Ecutoff=262 ± 34 keV. 

Considering NuSTAR together with INTEGRAL data 
for a temporally overlapping window we show that the 
broadband (3-150 keV) spectrum cannot be explained 
with a single spectral model. Extrapolating the NuS¬ 
TAR best-ht AGN model above 80 keV reveals a clear 
excess in INTEGRAL that we associate with emission 
from the jet. If we attempt to fit the combined INTE¬ 
GRAL and NuSTAR spectrum with a model including 
the AGN components and an additional power-law to 
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Energy (keV) 

Fig. 10.— Noise subtracted variance fraction of the mean for the 
four energy bands; 3-10, 10-20, 20-40, and 40—78keV. Statistically 
there is no evidence of a difference in variability as a function of 
energy. 



Fig. 11.— Light-curve showing GTI splits picked for periods 
where the flux remained roughly constant. 



Flux (10’'^ erg cm'^) 



Flux {10'’^ erg s'* cm'^) 



Fig. 12.— Flux in the respective bands 3-10, 10-20, 20-40, and 
40-78 keV as a function of a phenomenological pow er-law index 
fitted for each interval, GTI A-E, shown in Figure fTT] The unfilled 
points correspond to the interval GTI A. The solid line is a fit to 
all points, and the dashed line for intervals GTIB-E alone. 


describe the jet emission, we find strong degeneracies 
and no unique fit. We therefore fix the flux of the jet 
component from INTEGRAL 80-150 keV data and in¬ 
vestigate the range of fits to the broadband spectrum 
corresponding to the range of jet flux uncertainty. To 
further limit degeneracies we fix the relative reflection to 
R= 0.15. We find that the photon index for the coronal 
component ranges from 1.6 - 1.7, with a cutoff energy of 
30-70 keV. The jet photon index is poorly constrained, 
ranging from 0.5-1.5. 

Taking into account the high-energy y-ray spectrum 
from the Termi/LAT, and assuming a log-parabolic 
shape of the jet component, we find the local photon 
index at 100 keV of about 1.4. We caution that the com¬ 
plete lack of data in the 1-100 MeV range means that 
the broad-band shape of the jet component is uncertain. 
Nevertheless, the log-parabola model predicts a modest 
luminosity of the SED peak at 2 MeV. If the photon in¬ 
dex at 100 keV was in the range 0.5-1, the implied peak 
luminosity of the jet component would need to be signif¬ 
icantly higher, and that would be more demanding from 
the point of view of jet energetics. We note that a simi¬ 
lar log-parabola model a dded to an AGN comp onent was 
found to be successful by|Esposito et al.| (|2015|), who an¬ 
alyzed the Swift, INTEGRAL and Eermi data for 3C 273 
in different spectral states. In any case, we find that the 
jet component begins to dominate over the AGN compo¬ 
nent in the range from 30-50 keV. 


The cutoff energy range of 50-70 keV found with the 
jet included, is at the low end compared to other sources 


where a cutoff energy 

las been measured by NuSTAR 

(iBallantyne et al 

2014 

Brenneman et al. 20141 Marin- 

ucci et al.||2014 

datt et al. 

Balokoyic et al.pM I, 


to the poor constraints on the jet emission. If we exclude 
the jet component an upper limit for the energy cutoff is 
Ecutoff ~ 260 ± 40keV (see Table |^. 

We explore a physical model for the AGN component 
and fit with the coronal CompTT, which mod els spectra 


from a Gomptonizing coronal electron plasma (Titarchuk 
19941. If we fix the jet component parameters at the 


best-fit values given in Table we find an unusually low 
plasma temperature of kT = 12 ± IkeV, and high opti¬ 
cal depth Tp = 2.77 ± 0.06. Without the jet component 
included we find kT = 133^29 keV, and optical depth 
Tp = 0.33 ± 0.1. Since the jet contribution is poorly con¬ 
strained we can only conclude that the actual values lie 
between these extremes. 

The reflection fractions we derive for all models we 
consider are low (R= 0.02 — 0.2) (See Table |4 andj^. 
These indicate very weak reflection, as is oftennoundin 


Eracleous, Sambruna & Mushotzkyl |2000 

|/dziarski &: 

Grandi 

2001[ jSambruna et al. 2009 Bal 

antyne et al. 
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elude high inner disk ionizations 1 

Ballantyne, Ross & 

Fabian 2002 

), a change in inner disk geometry (|Era- 

cleous, Sam 

bruna & Mushotzkyj |2000l 

Lohhnk et al. 

2013), obscuration of the the central accretion how by 


the jet 1 

Sambruna et al. 

2009 

1 , 

black holes with retro- 

grade spin (Evans et al. 

2010 

), 

and dilution of the X- 

ray spectrum by jet emission ( 

Grandi, Urry & Maraschi 


20021. Since we still only find a small reflection traction 


even when taking into account the jet, it is likely not due 
to dilution of the jet. It has instead been suggested that 
a geometrically thick accretion di sk, which obscures part 


of itself from the external source (Paltani, Courvoisier & 


Walter 1998) could be responsible for the low reflection 


tractiorn Considering the possibility of a high optical 
thickness of the corona, it could also be that the corona 
itself is smearing out the reflection signature. 

3C 273 is well known for its variability. During 
the NuSTAR observation it went through several flux 
changes. Our study found no statistical evidence for a 
change in variability as a function of energy in the NuS¬ 
TAR band. The disentanglement and correlation of the 
jet and AGN component with respect to flux remains 
unclear. In the low energy band (3-10 keV) we find an 
inverse correlation between flux and slope, hardening the 
spectrum with increasing flux. This correlation seems to 
persist through all bands when excluding the highest flux 
bin, but appears to break down above 10 keV band when 
including the high flux bin. This inverse correlation in 
the low-energy band h as been previously ob served in the 
2-10 keV flux band by Kataoka et al. (2002) during their 
1999-2000 campaign with RATE. However, during their 
1996-97 campaign this correlation was no longer present. 
They attributed this change to the emergence of the AGN 
component during 1999-2000 when the source was at a 
higher flux than in 1996-97. We postulate that it may 
rather be that a higher flux level of the jet is responsible 
for the hardening of the spectrum. An inverse correla¬ 
tion is not common for AGN, where spectra typically get 
steeper with increasing flux (but we note that the oppo¬ 
site is generally true for blazars, both those associated 
with quasars, an d the lineless variety; r espective exam- 
ples are 3G279 ( iHayashi da et ah 2015), and Mkn 421 
(Balokovic et al.|2U13rr 'I'his could indicate that the 
variability is instead driven by the jet flux, hardening 
the spectrum as the jet flux contributes more strongly. 
The highest flux bin corresponding to GTIA does break 
the correlation, but considering the change in the light 
curve itself, there is an indication the source transitioned 
from one state to another at the end of GTI A. 

If we assume the two component scenario, it appears 
the AGN was bright at the beginning, clearly outshining 
the jet. Subsequently, lower flux softened the index at 
all energies, but it is not possible to tell what combina¬ 
tion of jet and AGN flux levels went into generating the 
variability sin ce the two compon ents are not expected to 
be correlated (|Soldi et al.||2008|). 

A long standing issue in 3G 273 is that the Comp- 
tonization of UV photons predicts a correlation to exist 
between UV and X-ray flux , which was claimed to have 
been observed in the past bylWalter h Gourvoisier (1992) 
but later found absent by ^ner^aE^'eF^ 1 21)07). In 
the jet/AGN component scenario the superposition of 
two uncorrelated model components would affect the re¬ 
liability of a correlation between UV and X-rays, and the 


disappearance of it may simply therefore be a matter of 
jet domination. 

The spectral structure inferred by the consideration 
of the NuSTAR and INTEGRAL data could possibly 
also be explained by invoking a two-component inverse 
Compton model arising in a relativistic jet. In such a 
scenario, the MeV-to-GeV range emission is likely to be 
dominated by scattering of the isotropic radiation exter¬ 
nal to the jet (as in External Radiation Co mpton, or 
ERG; see, e.g.|Sikora, Begelman & Rees|1994|), while the 
slight “hump” detected by NuS'IAR and modeled by us 
as Compton reflection from the accretion disk would be 
due to synchrotron self-Compton (SSC) emission, with 
the targe t photons being the synchrotron phot ons in the 
jet itself (Maraschi, Ghisellini & Celotti||1992 ). We con¬ 
sider such a scenario less likely: the NuSTAR data (as 
well as previous X-ray observations) detect a broad Fe K 
line, which is a signature of reprocessing from the accre¬ 
tion disk, which is likely to be accompanied by the reflec- 
tion component as in our modeling. Finally, as argued by 


Esposito et al. (2015), the timing analysis performed by 
those authors argues for the presence of the Seyfert com¬ 
ponent in addition to the jet component. However, any 
firm conclusions regarding the origin of the jet radiation 
- namely the relative dominance of the SSC vs. ERG 
component - are not possible with the data presented 
here. 


5. CONCLUSION 

In the coordinated observing campaign in 2012 we ob¬ 
served 3C 273 in a state where the jet flux was relatively 
weak. This made it possible to observe spectral sig¬ 
natures from coronal emission and constrain the weak 
Compton reflection in the hard X-ray band. We can sep¬ 
arate the spectral components from the inner accretion 
flow from the jet flux in the 3 - 150 keV band quantita¬ 
tively only by making an explicit assumption about jet 
flux level. To truly separate the two components and 
constrain the physical parameters of the corona, and/or 
address the possibility of SSC emission from the jet mask¬ 
ing as reflection, further simultaneous observations with 
NuSTAR, INTEGRAL and Eermi of the source in a state 
with low levels of flux from the jet are required. 

Eacilities: Fermi, CXO, INTEGRAL, NuSTAR, 

Suzaku, Swift, and XMM 
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